Interacting dark matter and modified holographic Ricci dark energy 
plus a noninteracting cosmic component 
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We investigate a spatially flat Friedmann-Robertson- Walker universe that has an interacting dark 
matter, a modified holographic Ricci dark energy (MHRDE), plus a third, decoupled component that 
behaves as a radiation term. We consider a nonlinear interaction in the dark component densities 
and their derivatives up to second order. We apply the % 2 method to the observational Hubble 
data for constraining the cosmological parameters and analyze the amount of dark energy in the 
radiation era for both MHRDE and holographic Ricci dark energy models. The former is consistent 
with the bound Q x (z ~ 1100) < 0.1 reported for the behavior of dark energy at early times while 
the latter does not fulfill it. 
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I. INTRODUCTION 

The current lore of modern cosmology indicates that 
95 % of the matter-energy in the Universe is harbored 
in the dark sector composed of dark matter (DM) and 
dark energy (DE) [H,U. Observations suggest that DM 
is a substantial nonbaryonic component that would be 
responsible for structure formation in the Universe Q. 
Likewise, DE has been the cornerstone of modern cos- 
mology ever since the measurements of supernova type 
la have confirmed its central role in explaining the cur- 
rent speeding up of the Universe [l[. The existence of 
DM and DE have been supported by observations such 
as cosmic microwave background and power spectrum of 
clustered matter among other probes [l|-|4[ . Here we will 
work within the framework of dynamical DE by taking 
into account as DE ansatz, p x , the modified holographic 
Ricci dark energy (MHRDE) based on the holographic 
principle [f|. The holographic DE model assumes that 
the quantum zero point energy density p\ = 3c 2 Mp/L 2 
is equal to DE density p x @, being L, an infrared cutoff 
that will be related with a cosmological length Q, Q, 

1, M, [? ]■ 

Despite the lack of microscopic theory to describe 
DM and DE, their existence offers some of the most 
compelling evidence for physics beyond the standard 
paradigm. However there is an unavoidable degeneracy 
between DM and DE within Einstein's gravity. There 
could be a hidden nongravitational coupling between 
them without violating current observational constraints, 
and thus, it is interesting to develop ways of testing dif- 
ferent kinds of interactions in the dark sector. 

The fraction of DE at recombination era should ful- 
fill the bound Q x (z ~ 1100) < 0.1 in order for the DE 



model be consistent with the big bang nucleosynthesis 
(BBN) data. Uncovering the nature of DE. as well as 
its properties to high redshift, could serve as an invalu- 
able guide to the physics behind the recent speed up of 
the Universe [12[. Current and future data were exam- 



ined for constraining the amount of dark energy at early 
times, and the data analyzed has led to an upper bound 
of S7 e < 0.043 with 95% confidence level (C.L.) in case 
of relativistic early dark energy, while for a quintessence 
type of EDE has given H e < 0.024, and although the EDE 
component is not preferred, it is also not excluded from 
the current data [12|- Another forecast for the amount 
of EDE was obtained with the Planck and CMBPol ex- 
periments [ill ; by studying the stability of the value 
il x (a ~ 10~ 3 ) ~ 0.03, it was found that la error coming 
from Planck experiment is a^ lanck ~ 0.004, whereas the 
CMBPol improved this bound by a factor of 4 (l3| . 

In the present work we investigate a Universe that has 
an interacting dark sector, and a decoupled component 
that could mimic a radiation term. We constrain the 
cosmic set of parameters by using the Hubble data and 
the severe bounds found for DE at early times. 



II. THE MODEL 

We will apply the holographic principle @ to a cosmo- 
logical model |6| associating the infrared cutoff L with 
the modified Ricci radius, thus we take L~ 2 in the form 
of a linear combination of H and H 2 . After that, the 
MHRDE p x = 3c 2 M 2 p L- 2 0, B, @, O, becomes 



a - (3 



H 



3a 



H- 



(1) 



a Elcctronic address: chimcnto@df .uba.arl 
b Electronic address: martin@df.uba.ar 



Here H = a/ a is the Hubble expansion rate, a is the scale 
factor and a,/3 are free constants. Introducing the vari- 
able r\ = ln(a/arj) 3 , with ao the present value of the scale 
factor and ' = d/di], the above MHRDE ([1]) looks like 



2 



a modified conservation equation (MCE) for both cold 
DM p c and the MHRDE with constant equations of state 
plus a continuity equation for the decoupled component 
named p m , 

p' = -ap c - PPx, p' m + {uj m + l)p m = 0, (2) 
after using the Fricdmann equation, 

P + Pm, 



3H 2 
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where p = p c +p x is total energy density of the interacting 



dark sector and p„ 
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-3(w m +i)_ j n connection with 



observations on the large scale structures, which seems 
to indicate that the Universe must have been dominated 
by nearly pressureless components, we assume that p c in- 
cludes all these components and has an equation of state 
(EOS) p c = while for p x and p m we have lu x = p x / p x 
and Lu m = Pm/pm, respectively. So that the whole con- 
servation equation (WCE) for the dark sector becomes 
p' = —p c — {lj x + l)p x . In addition, the compatibility 
between the MCE © and the WCE yields a linear de- 
pendence of the equation of state of the MHRDE with 
r = Pc/ p x , w x = (a — l)r + (/? — 1). Solving the linear 
algebraic system of equations © and p = p c + p x , we 
acquire both DE densities as functions of p and p' , 



Pp + p' 



p x 



ap + p 
A 



(4) 



with A = a — P, while the total pressure of the dark mix 
is p = Pc+Px = —p—p'- From now on we will use the sim- 
pler MCE (|2|) instead of the WCE with variable oj x , and 
introduce an interaction between both dark components 
through the term 3HQ into the MCE (J2J, 



p' c + ap c = -Q, p' x + f3p x 



(5) 



From Eqs. ((4]) and ([5]), we obtain the source equation 
for the energy density p, 



p" + (a + P)p' + afip = QA, 



(6) 



where the interaction Q between both dark components is 
nonlinear, and includes a set of terms which are the total 
energy density of the dark sector, its first and second 
derivatives (l4|. 
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Replacing Eq.© into Eq.© turns into a nonlinear 
second order differential equation for the total energy 
density p, a.k.a source equation (see [lH for further de- 
tails). Introducing the new variable y = p n into the 
latter equation one gets a second order linear differen- 
tial equation y" + (1 + 5)ny' = 0, whose solution is 

p = [Aa~ 3n(1+5 ^ +B] 1 ^ n , where A and B are positive 
constants. Note that Eq.© generalizes the EOS associ- 
ated with the modified Chaplygin gas in the sense that 
the exponents 1/n and — 3n(l + <5) arc different. Plugging 
p(a) in Eq. (U), we have both dark component densities 



A 



(7) 

being n, 5, and T positive constants. It turns out to 
be that the interaction (J7J gives rise to an effective fluid 
whose EOS includes all known Chaplygin gases [lH. In 
particular, it generates a new broad kind of gas that was 
called relaxed Chaplygin gas, [13], [13] ■ The aim of this 
work is to explore a model, which is characterized by an 
effective fluid having an EOS of a new kind of Chaplygin 
gas encoded in Eq. ([7]) by using observational constraints 
coming from the behavior of dark energy at early times. 
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while the total energy density is given by 



Pi 
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From Eq. pOl) we see that the Universe is dominated by 
radiation for uj m = 1/3 at early times where the dark 
components are negligible, after this epoch it enters an 
era dominated by dark matter and finally exhibits a de 
Sitter phase at late times. Then, the interaction allows 
a smooth interpolation between a cold matter era in the 
distant past (intermediate regime) to a speeding up stage 
at late times. Finally, the positivity of dark densities 
require to take P < 0, a > along with a > 1 + 6. From 
now on we adopt the latter restrictions. 



III. OBSERVATIONAL HUBBLE DATA 
CONSTRAINTS 

Using the observational H(z) data [13], we focus on 
the observational constraints on the parameter space of 
the MHRDE plus the third component. The statistical 
analysis is based on the x 2 function of the Hubble data 
which is constructed as (see e.g. [13]) 



■6)} (i-r) 
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where the 9 parameters are the cosmological set of pa- 
rameters, H b s (zk) is the observational H(z) data at the 
redshift Zk, and a(zt) is the corresponding la uncertainty 
[13]. Using Eqs. ©, ©, and © one gets the Hubble pa- 
rameter in term of the redshift x = 1 + z and the relevant 
cosmological parameters 
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FIG. 2: Plot of u) x , q, Q x , ^c, f^m in terms of the redshift z. 
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FIG. 1: Two-dimensional C.L. associated with la, 2a for dif- 
ferent 9 planes 



J[6] 



fl x0 A - a(l - Q. m0 ) 

(i + s)(i-n m0 ) 
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where we have used that 5(1 — J) = A J and the flat- 
ness condition: 1 = H x q + Q c q + Q m o- Now, we perform 
the statistical analysis by minimizing the x 2 function 
for obtaining the best-fit values of the random variables 
8 = {a, j3, S, n, fl x o, f2m0j -Ho}- The best-fit parameters 6 C 
are those values where Xmin(®c) leads to the local mini- 
mum of the x 2 (^) distribution. If the goodness condition 
(GC) xlo.f = X 2 mi n(0c)/(N - n) < 1 is satisfied then the 
data are consistent with the considered model H(z\8) 
We will estimate two parameters and the rest of 
them will be taken as priors, this will be done until all 
parameters have been varied [Fig.©], but for simplicity 
we will show some of these constraints only. As usual, 
in the case of two parameters, the 68.3%, 95.4% C.L. are 
made of the random data sets that satisfy the inequality 
Ax 2 = X 2 (d) ~ xLJfic) < 2.30, A x 2 < 6.17 [Fig.©]. ' 

The 2-dimensional CL obtained with standard \ 2 func- 
tion for two parameters is shown in Fig.© whereas 
the estimation of these parameters is briefly summa- 
rized in Table (|T|. In order to contrast the MHRDE 
ansatz with HRDE model, we performed the same es- 
timation in both cases. We find the best-fit at (n, Hq) = 
(0.64,73.60km s^Mpc" 1 ) with xl.o.f = °- 758 h Y usin S 
the same priors of the MHRDE ansatz except for the 
value taken by a. We also get the best-fit at (fl x o, Hq) = 
(0.76,73.81km s^Mpc" 1 ) along with xL./ = °- 760 
[cf. Table (JJ]. These findings show that both models 
are in agreement with the standard values of fl x o and 



H reported by the WMAP-7 project [T3|. Using the 
pair (n,H ) = (0.95, 74.09km s" 1 Mpc" 1 ) [Table ©], we 
show in Fig. © the behavior of the decelerating pa- 
rameter with the redshift; in particular, its present-day 
value is qo = —0.60, whereas for the HRDE model gives 
Qo = —0.59, so both values are in agreement with the ones 
reported by the WMAP-7 [l7|. In Fig.© we plot the 
density parameters f2 c , Q x , Q c , and find their present- 
day values for MHRDE model, namely, tt x o = 0.73, 
^mo = 10~ 8 , fl c o — 0.269, whereas HRDE ansatz leads 
to fl x0 = 0.79, fl m0 = 10~ 8 , and il c0 ~ 0.208. The 
EOS for dark energy does not cross the phantom barrier 
lo x = —1, so the model does not exhibit a quintom phase 

m 

We devote our attention to the behavior of DE at 
early times within the framework of interacting dark 
sector plus a decoupled third component. We begin by 
setting a = 4/3 to get a radiation contribution in the 
total density, and, as is well known in the radiation 
era the stringent bound VL x (z ~ 1100) < 0.1 should be 
fulfilled by the model in order to be consistent with 
data coming from the recombination era. In the case 
of the HRDE model with (a, f3, S, n, £I x q, fi m o, Hq) = 
(4/3, -0.01, 0.01, 0.8, 0.76, 10~ 8 , 73.81 km s^Mpc" 1 ) 
the fraction of DE is tt x (z ~ 1100) = 0.24; there- 
fore, it does not satisfy this bound. For the 
MHRDE case, we take (a, (3, S, n, fl x o, fi m o, H ) = 
(1.01, -0.01, 0.01, 0.9, 0.71, 10" 8 , 73.94 km s" 1 Mpc" 1 ) 
along with x\ o / = 0.763 < 1. The latter case leads 
to an EDE with Q x (z ~ 1100) = 1.2 x 10~ 8 < 0.1 
that is consistent with both the bounds reported in 
(Tj , and with the future constraints achievable by 
Planck and CMBPol experiments [l~3| because at large 
rcdshifts we get some bounds for the DE such as 
fl x (z ~ 10 10 ) = 1.3 x 10~ 27 , fl x (z ~ 10 16 ) = 5.8 x 10~ 44 , 
and fl x (z ~ 10 20 ) = 7.23 x 10 ~ 55 . Then, taking into 
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account the third component as a radiation term or 
a nearly radiation contribution helped to validate our 
model concerning the present-day constraints, but also 
it shows that the value of the cosmological parameters 
selected are consistent with BBN constraints. 



IV. CONCLUSION 

We have considered a model of three components com- 
posed of an interacting dark sector with MHRDE and a 
nonintcracting component. We have introduced a non- 
linear interaction such that the dark mix interpolates be- 
tween a cold DM regime at early times (after a radiation 
era) and a de Sitter phase in the far future. 

In the case of 2D C.L., we have found that the 
MHRDE model leads to the best- fit at (Q x0 ,H ) = 
(0.71,73.94km s^Mpc" 1 ) with xlo.f = °- 763 [ Fi S-© 
and Table fT}] wni l e tnc HRDE leads to a best-fit 
at (O x0 ,i? ) = (0.76, 73.81km s" 1 Mpc" 1 ) plus a 
Xd o f — 0.760 < 1. Both models arc in agreement 
with the standard values of tt x o and Hq reported by 
the WMAP-7 project [13 and fulfill the GC. Using 
the best fit-value (n,H ) = (0.95, 74.09km s" 1 Mpc" 1 ) 
[ Table (HJ], we have shown that the behavior of 
the q{z) exhibits an accelerated phase at late times 



[see Fig. (|5J)] and obtained that its present-day 
value is = —0.60, whereas for the HRDE model 
yielded qo = —0.59, so both values are i n ag reement 



with the ones reported by the WMAP-7 [17J- Since 
the HRDE model with (a, (3, S, n, £l x o, ^ m o, Hq) — 
(4/3, -0.01, 0.01, 0.8, 0.76, 10" 8 , 73.81 km s" 1 Mpc" 1 ) 
along with xlof = °.760 < 1, leads to an EDE 
Q x (z ~ 1100) = 0.24, then this model docs 
not fulfill the bound for EDE at the recombina- 
tion era. On the other hand, for the MHRDE 
case we have taken (a, f3, <5, n, Q x q, fimd Hq) — 



(1.01,-0.01,0.01,0.9,0.71,10- 
along with Xd.o t = 0.763 < 



,73.94 kms^Mpc' 1 ) 
D y — u.iuj ^ 1. The latter case leads 
to an EDE Q x (z ~ 1100) = 1.2 x 10~ 8 < 0.1 which is 
consistent with both the bounds reported in fl2j . and 
with the future constraints achievable by Planck and 
CMBPol experiments [l3j . 
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J 



2D Confidence Level 


Priors 


Best fits 


Xdof 


(a, /3,S, f) x0 ,fimo) = (1.01, -10- 2 ,10- 3 , 0.73, 10" 8 ) 


(n,H ) = (0.95,74.09) 


0.766 


(a,/3,5,n,n m0 ) = (1.01, -lO" 2 , 10~ 3 , 0.9, 10" 8 ) 


(Q x o,Ho) = (0.71,73.94) 


0.763 


(a, P, 5,n x0 ,H ) = (1.2,-10- 2 ,0.1,0.73,74.20) 


(n,Sl m0 ) = (0.63, 10" 8 ) 


0.764 


(a, 8, n, n x0 , Q m0 ) = (1.26, 0.1, 0.8, 0.73, 10" 8 ) 


(p,H ) = (-0.09,74.09) 


0.771 


(a, p, Q x0 , Q mQ , Ho) = (1.2, -10~ 2 , KT 3 , 0.73, KT 8 , 74.2) 


(S,n) = (0.1,0.63) 


0.764 


(/3, n, Q x0 , n m0 , H ) = (-10- 2 , 0.4, 0.73, 10" 8 , 74.2) 


(6, a) = (0.25, 1.37) 


0.766 



TABLE I: Observational bounds for the 2D C.L. obtained in Fig. (JTJ by varying two-parameters. 



